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We analyze the influence of unconventional superconductivity on the magnetic excitations in 
the heavy fermion compound UPd2Al3. We show that it leads to the formation of a bound state 
at energies well below 2Ao at the antiferromagnetic wave vector Q=(0, 0, 7r/c). Its signature is 
a resonance peak in the spectrum of magnetic excitations in good agreement with results from 
inelastic neutron scattering. Furthermore we investigate the influence of antiferromagnetic order 
on the formation of the resonance peak. We find that its intensity is enhanced due to intraband 
transitions induced by the reconstruction of Fermi surface sheets. We determine the dispersion of the 
resonance peak near Q and show that it is dominated by the magnetic exciton dispersion associated 
with local moments. We demonstrate by a microscopic calculation that UPd2Al3 is another example 
in which the unconventional nature of the superconducting order parameter can be probed by means 
of inelastic neutron scattering and determined unambiguously. 



I. INTRODUCTION 



The relationship between unconventional superconduc- 
tivity and magnetism in heavy-fermion systems and tran- 
sition metal oxides is one of the most interesting research 
areas in condensed matter physics. In both cases it is 
widely believed that the magnetic degrees of freedom 
play an essential role in the formation of superconduc- 
tivity. Furthermore, unconventional superconductivity 
yields strong feedback on the magnetic spin excitations 
in these systems below the superconducting transition 
temperature T c . One example is the famous so-called 
resonance peak observed in high-T c cuprates by means of 
inelastic neutron scattering (INS) 1 whose nature is still 
actively debated^ Remarkably, it has been found that 
INS reveals the formation of a new magnetic mode in 
the superconducting state of the uranium based heavy- 
fermion compound UPd2Al3 with T c = 1.8 K. 3 Its sharply 
peaked intensity, its temperature dependence and the en- 
ergy position well below 2Ao (with Ao being the maxi- 
mum of the superconducting gap) strongly resembles the 
resonance peak seen in high-T c cuprates. This is par- 
ticularly remarkable, since the origin of superconduc- 
tivity in cuprates and UPd2Al3 seems to be different. 
While frequently discussed scenarios in cuprates are a 
spin-fluctuation mediated Cooper pairing or the electron- 
phonon interaction, in LT^Ahj a magnetic-exciton me- 
diated pairing model has been proposed 5 - based on avail- 
able experiments. The latter model is built on the dual 
nature of the 5/ electrons. It consists of localized 5/ 2 
crystalline electric field (CEF) states which disperse into 
a magnetic exciton band due to intersite interactions and 
a conduction electron band* formed by itinerant 5/ elec- 
trons with enhanced hybridization. The model success- 
fully explains the formation of unconventional supercon- 
ductivity in this compound 5 based on the virtual ex- 
change of the magnetic excitons between itinerant quasi- 
particles. 



It is important to note that the resonant spin excita- 
tions in superconducting cuprates can be seen as a direct 
consequence of the d x 2_ y 2-wsive symmetry of the super- 
conducting order parameter. Namely, the resonance peak 
occurs only if the order parameter changes sign in the 
first Brillouin zone (BZ)£ Thus, INS can be considered 
as a bulk probe for the unconventional nature of super- 
conductivity in these compounds. Therefore it is impor- 
tant to search for such an effect in other unconventional 
superconductors as well. In this paper we analyze the 
consequences of the unconventional pairing on the mag- 
netic excitations in UPd2Al3. We show that in addition 
to the magnetic exciton dispersion present in the nor- 
mal state, unconventional superconductivity induces the 
formation of a bound state below T c with an associated 
resonance peak in the magnetic spectrum at the antifer- 
romagnetic (AF) wave vector Q=(0, 0,7r/c) where c is a 
lattice constant along the crystallographic z axis. Its fre- 
quency is well below 2Aq and in good agreement with 
experimental data. We show that similar to cuprates 
the resonance peak in UPd2Al3 is a consequence of an 
unconventional superconducting order parameter which 
changes sign at regions of the Fermi surface connected 
by the antiferromagnetic wave vector Q. We analyze the 
influence of antiferromagnetism on the formation of the 
resonance peak and surprisingly find that its intensity is 
enhanced due to the reconstruction of the Fermi surface. 
We find that the dispersion of the resonance peak away 
from Q is controlled by the momentum dependence of 
excitations of the localized magnetic moment (magnetic 
exciton) . 



The resonance peak in UPd2Al3 has also been stud- 
ied by Bernhoeft et alZ within a phenomenological two 
component spin susceptibility model. However, in our 
microscopic calculations we show that antiferromagnetic 
order plays a crucial role in the formation of the reso- 
nance peak below T c . 
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II. THE HAMILTONIAN 

Following previous consideration by McHale et al.— we 
use the low-energy Hamiltonian describing the interac- 
tion of the itinerant / electrons and magnetic excitons 
originating from localized 5/ 2 crystalline electric field 
(CEF) states: 

H o = X^p c p^ c p ct + ( a q a q + 2) 

per Q 

- jf^ c \> a v z a i3 c P+<i,pK (a n + at q ) , (1) 
p.q 

where A 2 = A ^ EF , and Aqef = 6 meV is the energy 

gap of the 5/ 2 electrons between the ground and first 
excited states in the crystalline electric field. The disper- 
sion of the magnetic excitons is approximately described 
by uj(q z ) — io ex [\ + (3 cos(cg z )] with < /3 ~ 1, where g 
is the coupling constant between the itinerant electrons 
and the localized magnetic moments. We adopt param- 
eter values (jjp X = 5.5 meV and (3 = 0.72. Note, here we 
follow Ref. in assuming that only the a z component 
of the conduction electrons can excite magnetic excitons. 
Therefore the spin-space isotropy is broken in a maxi- 
mal (Ising) way. As a result the usual classification of 
Cooper pairs into spin-triplet and spin-singlet states is 
not valid and the notation equal and opposite spin pair- 
ing states should be better used instead. However, we will 
still speak of singlet and triplet Cooper-pairing states as 



commonly done. 

Eq. (JlJ gives rise to fermionic and bosonic self-energies 
and is particularly relevant for electron-hole states sepa- 
rated by the antiferromagnetic wave vector Q=(0,0,7r/c). 
Previously it has been shown that this interaction ex- 
plains superconductivity in UPd2Ai3j£ We define the 
electron and magnetic exciton Green's functions as fol- 
lows: 

G CT<T ,(p,iu, m ) = -(T T c pCT (T) C t (0)) , 

\ / FX 

D(q,iv n ) = - (T TOq (r)oL q (0)) , (2) 

\ / FT 

where a q (r) = a q (r) + aL (t) and D is essentially the 
pseudospin susceptibility. The bare magnetic exciton 
Green's function is given by D a (q, iv n ) = — A( ^ EF ^ l-^ , 

Due to the interaction of the magnetic excitons with con- 
ducting electrons, the feedback effect on the former re- 
sults in 

D = -Jh = ^ ( 3 ) 

1 - D n u 2 - uj\ + 2 Wq n ' 1 ' 

where the magnetic exciton self-energy is given by 

tt / - \ 2^CEF / ■ \ ,.s 

n Q {q,iuj n )=g xo(q,«w„) . (4) 

Here, the spin susceptibility of the conduction electrons 
in the superconducting state is 



Xo(q,iw„) = -- ^2 [ G ( k ' w m )G t (k + q, ioj m + iu n ) + F(k, iuj m )F^(k + q, iuj m + iuj n )] , (5) 

where bare Green's functions of superconducting electrons are G(k, iu m ) = rrj&rki » f(k,iu m ) = , Jtf, a-i ■ A 

straightforward evaluation of the sum over the Matsubara frequencies gives (at T = K) 

Im *°<"-"> = w / dS " ( l - &±3 |^ir t±3 ) '<» - E ^ - E » ■ < 6 > 

I 



The Fermi surface of the conducting electrons is almost 
like a cylinder with weak dispersion along the z direction. 
Neglecting the anisotropy of dispersion in the plane yields 
Ck = £fei + e kz - fi = e ± w 2 + 6|j cos(cfcz) - /u (e\\ < e ± , 
w = kj_/k < 1) and Ek = yfk + A 2 .. Here, we ap- 
proximate the hexagonal unit cell by a circle with radius 
ko chosen so that the hexagon and the circle have the 
same area. Furthermore, we assume a parabolic disper- 
sion in the plane. Due to the Ising-type anisotropy of the 
interaction between conduction electrons and magnetic 
excitons it has been previously founds that both pure 
paramagnetic, i.e., spin-singlet states -4= (| Tl) — | IT}) 



with cos(ck z ) momentum dependence and spin-triplet 
states (S z = 0) ^= (| Tl) + I IT)) with sm(ck z ) have the 
highest (degenerate) superconducting transition temper- 
ature. Thus, in the following we will consider the two 
superconducting order parameters, A£. = Ao cos(ck z ) 
and A^ = Ao siri(ck z ) as the most relevant ones in this 
model. 

Let us now discuss the consequences of the behavior of 
Irn^o f° r the magnetic exciton dispersion which follows 
from Eq. ([3]) . The dispersion of the magnetic exciton in 
the presence of coupling to the conduction electrons is 
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given by 

u? = uj\ - 2 5 2 A C E F Re Xo(q, u) . (7) 

In the normal paramagnetic state, Rexo is a constant 
at low frequencies determined in our case by the curva- 
ture of the Fermi surface along the k z direction. At the 
same time Im^o — *7^: where 7 is a Landau damp- 
ing constant. Thus, the bare magnetic exciton acquires 
a linewidth and renormalizcs slightly by a certain con- 
stant which changes the original position of cj q down- 
wards. In the superconducting state the renormalization 
is strongly dependent on the gap symmetry. For example, 
in the conventional s-wave state Re xo is less than its nor- 
mal state value. First, the s-wave superconducting gap 
gives a negative contribution to xo as follows from the 
coherence factor in Eq. ©, and second it yields the spin 
excitation gap structure in lmxo as follows from the S 
function. Therefore, the feedback of the conduction elec- 
trons on the magnetic exciton becomes weaker and yields 
a shift of the magnetic exciton dispersion towards higher 
frequencies in the superconducting state with respect to 
its normal state value. 

At the same time, the unconventional character of the 
superconducting gap entails immediate consequences for 
Eqs. © and ([7]). Namely, at the Fermi surface (i.e., for 
Ck = 0) one finds A* +Q = -A£ as well as A£" +Q = -A^ r 
where Q is the antiferromagnetic wave vector. Thus, 
in both cases the anomalous coherence factor equals 
2 for all k z momenta which is in strong contrast to 
the usual s-wave symmetry of the superconducting 
gap. Simultaneously, the 5 function in Eq. ([6]) starts 
to contribute for £l cr (k z ,q z ) = Min (-Ek+zg, + -E-k) = 

«/ejj [cos(cfc z + cq z ) -cos(ck z )] 2 + (|A fe J + \A kz+qz \) 2 

which is determined by the value of the gap, Ak and 
the quasiparticle dispersion along k z axis, en. Note 

that £l s cr (k z ,q z = ir/c) = 2 1 cos(cfc z ) | ^/Aq + ejj and 

Sl%{k z ,q z = n/c) = 2^ + (e| - Ag) cos 2 (cfc 2 ) for 

the singlet and triplet superconducting gap, respectively. 
It is interesting to note that for the singlet case there 
is no gap in Im^o an d therefore, it is nearly the same 
as in the normal stated [see Fig. [TJa)]. Then Kexo is 
a constant at low frequencies. Therefore the magnetic 
exciton dispersion will be simply shifted in proportion to 
the total value of Re%o exactly as in the normal state. 
Correspondingly, no strong feedback on the magnetic 
exciton due to superconductivity takes place. 

However, one sees that in the case of the triplet order 
parameter Im xo will be gapped at least up to values of 
fi*^ = |2Ao|(eii ^ Ao). Then due to a combined effect 
of the anomalous coherence factor and the S function, a 
discontinuous jump in Im^o occurs at about u> = 2Ao for 
the triplet order parameter. Via Kramers-Kronig trans- 
formation the discontinuous jump in Im^o yields a log- 
arithmic singularity in Re^o- Note, the logarithmic sin- 
gularity in Fig. l.(b) is suppressed by a weak damping. 
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FIG. 1: (Color online) The (a) imaginary and (b) real parts of 
the longitudinal component of the conducting electrons spin 
susceptibility xo(Q,<^) as a function of frequency for singlet 
(red dashed curve) and triplet (black solid curve) Cooper pair- 
ing. The curves for the normal state (not shown) almost co- 
incide with those for the spin singlet superconducting state. 
Here and in the following we use T = 0.15 K, Ao = 0.5 meV, 
e 11 = 25 meV and a damping parameter 35 fieV for the nu- 
merical calculations. 



Furthermore, below il*' r the Re xo is increased with re- 
spect to its normal state value for w ^ as shown in Fig. 

raw- 

A frequency dependence of the HexoiQ^) can yield 
more than one solution of Eq. (7). In order to demon- 
strate how those solutions can be found above and below 
T c we illustrate in Fig. 2 the possible characteristic be- 
havior of cj 2 — 2g 2 AcEF Rexo(q, w). Here, we assume 
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FIG. 2: (Color online) Illustration of the solutions of Eq. (jTj) 
at wave vector q = Q. In the normal state there is only one 
crossing point between the lo 1 curve (solid) with the w q line 
(dotted) yielding the frequency of the magnetic exciton. Note 
that u q may slightly differ from the bare exciton dispersion 
due to Rexo(Q,<^) = const in the normal state. In the su- 
perconducting state due to the strong frequency dependence 
of Rexo(q,w) (and/or Rello) one finds several intersecting 
points of ui 2 with cj^ — 2g 2 AcEFRe Xo(Qj w )- The lowest pole 
[u> T < 2Ao) occurs at very small damping (Im^o is zero or 
small) resulting in a resonancelike peak in ImZ)(q,a;). The 
second crossing point is not visible in Im D due to a large peak 
in lmxo or strong damping around 2Ao. The third crossing 
point, uJm occurs at energies larger than 2Ao and represents 
the feedback effect of superconductivity on the magnetic ex- 
citon. 
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that in normal state the Re%o is almost frequency in- 
dependent and magnetic exciton's peak position shifts 
slightly in superconducting state with respect to its nor- 
mal state value. Due to the gap structure in supercon- 
ducting state, depending on the value of g, a new pole 
may occur at energies less than 2Ao. If Im%o is small or 
zero at these frequencies, the total Imfl, i.e., the spec- 
tral function of magnetic excitations, shows a resonance 
peak which occurs only in the superconducting state. 
This agrees well with experimental INS data 3 . Moreover, 
at higher energies one finds in addition two more poles 
in Eq. ([7]). The latter yields an additional structure in 
ImD which is a renormalized magnetic exciton with fi- 
nite damping. This typical behavior of the susceptibility 
can be found in Fig. 3. 

So far we have ignored the coexistence of antiferro- 
magnetism and superconductivity in UPd2Al3 for the 
conduction electrons. 8 Antiferromagnetic order results in 
UPd2Al3 due to the interaction of neighboring uranium 
ions. This leads to a dispersion for exciton state and 
eventually to an antiferromagnetic instability and a new 
ground state.— The unit cell is doubled and the Brillouinc 
zone is correspondingly reduced. The new dispersion of 
the conducting electrons enters the expression for the 
spin susceptibility. Then the solutions of Eq. (7) must 
be redetermined. This is done in the following. 

The total Hamiltonian is 

H = Ho + m^crc^ +Q(T Cp ff (8) 
p 



where to denotes the value of the effective antiferro- 
magnetic staggered field. This term leads to a splitting 
of the quasiparticle energy dispersion into two bands.— 
In particular, the Hamiltonian ([5]) can be easily diag- 
onalized by a unitary transformation^ and the result- 
ing energy dispersions are E^ = J (e k ) 2 + A 2 , with 

= [e^±yj (e^) 2 + to 2 ]. Here, we have introduced = 

3 (6k + 6k+q) and e k = 1 (& - Ck+q). The AF Fermi 
surface of the two bands consists of two disjoint cylin- 
ders in the reduced AF BZ \p z \ < \ (see Ref. H3). In a 
pure AF state, the real part of the susceptibility xo at low 
energies is determined by the intraband processes and can 

be approximated by Xo(q,^) °c ^k 4,5^7^ 

where A k „ is the AF coherence factor. At wave vec- 
tor Q for u> = the susceptibility is proportional to the 
density of states and decreases rapidly to zero as one in- 
creases frequency. This is a consequence of the equality 
= £ k+ q- Therefore the renormalization of the mag- 
netic excitons due to conduction electrons can be safely 
ignored. 

Most importantly, in the superconducting state coex- 
isting with AF, the imaginary part of the spin suscep- 
tibility of the conduction electrons including intraband 
and interband scattering is given for T=0 K, q = Q, and 
lu > by 



Imxo(Q,w) 



k 



A k 



2m 2 



£z%-i£-i?k) i 



e e 



k c k 



"A k 



2m 2 



2 < 



(9) 



Here the first and the third terms describe the intraband 
quasiparticle pair creation while the second term refers to 
the corresponding interband process. Note that Eq. © 
contains two types of coherence factors, i.e., due to super- 
conductivity and antiferromagnetic order, respectively. 
As usual, the low-energy behavior of Imxo is dominated 
by the intraband contributions. We also assume that 
the presence of antiferromagnetism does not change the 
qualitative behavior of the superconducting gap, i.e., the 
position of the line node and the corresponding change of 
sign of the superconducting order parameter remain the 
same although some higher harmonics may appear ^ As 
in Eq. ([6]), the superconducting coherence factors equal 2 
for k z momenta close to the Fermi surface. At the same 
time, the reconstructed conduction bands in the AF state 



have only one-half the original period, i.e., = £ k+ Q- 
Therefore all parts of the Fermi surface can be connected 
by the antiferromagnetic wave vector Q and simultane- 
ously have a sign change A k+ Q = — A k of the gap except 
at the nodal points. As a result, Imxo(Q, w ) is non-zero 
for uj > due to the contribution of the nodal states 
both in the singlet and the triplet Cooper-pairing cases. 
With increasing frequency Im xo increases up to energies 
of about 2A and then decreases [see Fig. [3ja)]. The 
functional dependence of Im xo at low frequencies resem- 
bles the behavior of the density of states except that the 
structure occurs at around 2Ao. Correspondingly, the 
real part of xo(Q, w = 0) is the same as in a pure AF 
state. However, away from u) = it does not drop as in 
the pure AF state but increases quadratically up to about 
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FIG. 3: (Color online) Results for Xo(Q,<^) and D(Q,w) 
in the coexisting AF and singlet superconducting state, 
(a) Calculated real and imaginary parts of the longitudi- 
nal component of the conduction electron spin susceptibil- 
ity, xo(Q,u)- (b) The real and imaginary part of the 
total pseudospin susceptibility D(Q,u>). Here m = 50 
meV, g — 10 meV. We use the gap function Ak tx 
Ao [cos(cfc z ) — | cos(3cfc z ) + -jj cos(5cfc z )] ( Ref. ). 

2A due to the structure of Imxo induced by the super- 
conducting gap. Only then does Re;\;o drop to small 
values. Altogether Rexo(Q, w ) increases in the super- 
conducting state for uj < 2Ao due to the unconventional 
nature of the superconducting order parameter. How- 
ever, the pure resonance (bound state) in ImD is not 
realized due to finite damping. An additional pole in 
ImD still exists at frequencies smaller than 2Aq due to 
a strong increase of Re \o &t small frequencies as shown 
in Fig.^b). At higher frequencies it becomes small and 
thus another pole appears corresponding to the broad- 
ened original magnetic exciton. Thus, ImD has a two- 
pole structure as shown in Fig. [3fb). Note, in order to 
increase the intensity of the low-energy pole in Im D we 
include the contribution of the higher harmonics to the 
gap function. As already mentioned, they are due to the 
presence of AF order—. 

Finally we discuss the dispersion of the magnetic exci- 
tations away from Q along the q z direction. In Fig. [4] we 
show the calculated momentum and frequency dependen- 
cies of ImD (q z , cu). It is clear that as soon as q z ^ir/cthe 
original magnetic exciton has a strong upward dispersion 
in the normal state. Therefore, effects connected with 
renormalization induced by superconductivity will also 
be shifted towards higher energies. The pole induced by 
superconductivity shows dispersion similar to the mag- 
netic exciton [see Fig. 0]. For both singlet and triplet 
order parameters our results are in fair agreement with 
recent INS datap 1 ^ Namely, in the superconducting state 
one finds two distinct energy dispersions, one being the 
resonancelike feature with high intensity inside the super- 
conducing gap and the second, that of localized magnetic 
excitons renormalized by the conducting electrons. An- 
other interesting point worth noting is that due to the 
doubling of the unit cell and the equality = £^ + q, 
the effect of the sin(cfc z ) and cos(cfc z ) gaps leads to a 
very similar behavior for Imxo- The slight difference in 




o) (meV) 



FIG. 4: (Color online) Contour plot of the imaginary part of 
the total pseudospin susceptibility as function of frequency, ui 
and q z momentum. One clearly observes two distinct peaks 
at Q. The one at low energies represents the resonance peak 
induced by the feedback of superconductivity and the one at 
higher uj is the renormalized magnetic exciton. Away from Q 
both peaks disperse upward in energy following the behavior 
of the normal state magnetic exciton. 

the absolute magnitude arises from the different densi- 
ties of states at those regions of the Fermi surface where 
the maximum of the singlet and the triplet gaps occur. 
Altogether this does not change the functional form of 
Imxo- We mention that thermal conductivity results 
in a rotating field 14 are compatible with both cos(ck z ) 
and sin(cfc z ) order parameters while the observed Knight 
shift 15 seems to favor the former. Interestingly, we also 
found that a recently proposed superconducting gap with 
cos(2cfc 2 ) symmetry 1 ^ does not lead to the formation of 
low-energy spin excitations around wave vector Q in the 
superconducting gap. The reason is that its momentum 
dependence yields no change of the sign of the super- 
conducting order parameter, Ak = Ak+Q, and thus no 
constructive contribution can result from the anomalous 
coherence factor. 

In conclusion, we have investigated the effects of super- 
conductivity on the magnetic excitations in the uncon- 
ventional superconductor UPd2Al3. In particular, due to 
the change in sign of the superconducting order parame- 
ter the conduction electron susceptibility is enhanced in 
the superconducting state which yields an additional pole 
(bound state) in the total susceptibility. We further an- 
alyzed the role played by antiferromagnetism and found 
that its presence increases the spectral weight of the res- 
onance due to the doubling of the unit cell. However, 
the resonance peak in the AF phase becomes a virtual 
bound state due to finite damping. Finally we point out 
that UPd2Al3 is another known example where the un- 
conventional nature of the superconducting order param- 
eter yields a structure in the magnetic susceptibility as 
in layered high-T c cuprates. Therefore it can be regarded 
as a model system of unconventional superconductivity 
studied by inelastic neutron scattering. 

The authors acknowledge helpful discussions with G. 
Zwicknagl, J. -P. Ismer, and A. Klopper. 
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